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Abstract 

Taking advantage of the conformal equivalence of f{R) theories of gravity with 
General Relativity coupled to a scalar field we generalize the Israel junction conditions 
for this class of theories by direct integration of the field equations. We suggest a 
specific non-minimal coupling of matter to gravity which opens the possibility of a new 
class of braneworld scenarios. 
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typeset using pTpT^.ds (Ver.0.9) 



§1. Introduction 



In 1918 Weyl^) was the first to consider Lagrangians for gravity wliicli are a linear com- 
bination of the three quadratic scalars i?^, RabR'^^ and RabcdR^^^'^ formed out of the 
scalar curvature, the Ricci and Riemann tensors of a metric qablI- 

Pauli^-* and Eddington^-' then noticed that the Schwarzschild metric was a solution of 
the corresponding vacuum field equations. Since at that time all precision tests of General 
Relativity theories relied on the Schwarzschild metric, these authors concluded that quadratic 
theories were a priori just as viable as ordinary General Relativity based on the Einstein- 
Hilbert Lagrangian R. (As noted by Chiba et al.J^ this seems to be still a fairly widespread 
belief.) 

The question of the unicity of the Schwarzschild metric in such theories in vacuo was first 
investigated by Buchdahl^^ who showed that the Schwarzschild metric was the only spher- 
ically symmetric, asymptotically flat, solution of the vacuum pure R^ field equations. The 
subject was considered anew much later by Whitt^-* (in the case of the R + a2R^ Lagrangian 
in 4 dimensions), and Mignemi and Wiltshire^"-' who showed that the Schwarzschild metric 
was the only D-dimensional, spherically symmetric, asymptotically fiat, vacuum solution 



with a regular horizon, for all polynomial f{R) = R + Yl o-nR^ with 02 > Q**-'. 

Now, since the f{R) field equations are fourth order differential equations for the metric, 
they possess "runaway" solutions on top of solutions which smoothly tend to solutions of the 



Einstein equations in the limit f{R) H***\ The question of whether the Schwarzschild 



metric is the f{R) solution outside a distribution of ordinary matter (rather than a black 
hole), either point-like or extended, must therefore be raised. Pechlaner and Sexl^^^ showed 
that, in fact, in pure R^ theory the metric cannot be asymptotically fiat as soon as the field 
equations have a right-hand side describing matter with positive energy density. They also 
showed that in R + a2-R^ theory the metric can be asymptotically fiat but that, at linear 
order around Minkowski spacetime, it is not the linearized Schwarzschild metric. The origin 
of such results is clearly explained by Havas^^-* : The Green function for the (second order) 
Einstein equations, which at lowest order reduce to V^Gat = — 5(r), takes the familiar form 
Gat = which yields the linearized Schwarzschild solution. On the other hand, in the case 
of the (fourth order) pure R^ theory the Green function solving the linearized field equations 



*) The fact that the RabcdR^^^'^ can be traded for — (i?^ — ARabR^^) (in 4 dimensions) was appar- 
ently first shown by Bach, and then by Lanczos^' (the Gauss-Bonnet theorem'^)). For a historical review, 
see 4). 

**•* As emphasized in 10), other solutions exist if one relaxes the condition of asymptotic flatness and 
allows for asymptotically de Sitter or anti-de Sitter spacetimes. 
***) rpj^-g ^g^g used by Starobinsky^^^ to build the first infiationary cosmological model. 
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V^G = —S{r) is G = r/2, which yields a divergent metric. As for Gn = 1/r it satisfies 
V'^Gn = — V^5(r) ; this means that in pure theory the source for the Schwarzschild 
metric is not a delta function but its second derivative which does not represent a point-like 
distribution of matter with positive energy density. Finally the Green function for R+a^R^ / 6 
theory is G = (1 — e^'"/'*)/r, yielding the Pechlaner-Sexl metric (after correcting a couple of 
typos in their equation (23)). (See also 6), 14).) 

This question of unicity or non-unicity of the Schwarzschild (or Schwarzschild-(A)dS) 
solution was revived recently when f{R) theories were invoked in an attempt to explain the 
observed present acceleration of the universe by means other than a cosmological constant ^^-^ 
(see also 16) and references therein). In that context the question became whether or not 
the Schwarzschild (or Schwarzschild-de Sitter) metric is, at least approximately, a solution 
of the /(/?) field equations in the presence of localized sources such as the Sun. In 17), 
static spherically symmetric solutions for the special case /(/?) = R — fi^/R were built 
numerically with matter represented by a perfect fiuid (see also references in 17)). It was 
found that, if the metric tends to an appropriate de Sitter limit to explain the acceleration 
of the universe then the PPN parameter 7 measuring, e.g., the light bending by the Sun is 
of order 1/2 instead of 1,-*^^) which rules out these "dark energy" models. In 7) and 19) (see 
also references therein), the field equations were solved at the linear approximation around 
Minkowski spacetime and similar results were found for a wide class of f{R) theories. 

In many of the above mentioned analyses advantage is taken from the fact that f{R) field 
equations are conformally equivalent to Einstein gravity with a minimally coupled scalar field 
(this was first shown by Higgs'^'^' for f{R) = R^ and by Teyssandier and Tourrenc^^^ in the 
general case). This property puts the theories on a more familiar footing but, in itself, does 
not modify either the mathematics or the physics of the problem. Nevertheless conclusions 
concerning what are the "correct" solutions of the field equations are still controversial. In 
particular the validity of the linear approximation has been challenged (see e.g. 22) and 
references therein) and the question of how one recovers the (de Sitter-) Schwarzschild 
solution in the Einstein hmit does not seem to be settled yet (see e.g. 23) and references 
therein) . 

In this paper, taking advantage of this conformal equivalence, we investigate junction 
conditions for a brane (i.e. an infinitesimally thin domain wall) in f{R) theories. The topic 
has already been investigated in 24), where, however, the questions of the Einstein limit and 
the coupling to matter on the brane were not addressed. Here, as shown below, we are able 
to see clearly the possible irregularity which may appear in the Einstein limit f{R) —>■ R- 
Then we suggest a specific non-minimal coupling of matter on the brane to gravity which 
opens the possibility of a new class of braneworld scenarios. 



3 



The paper is organized as follows. In Sec. [21 we consider general quadratic theories of 
gravity and discuss possible freedom in the choice of the junction conditions depending on 
how singular one allows the metric to be. Then, as an example, we formulate the junction 
conditions by requiring the metric to be least singular, namely its first and second derivatives 
are continuous across the brane. Then, in Sec. [3l focusing of f{R) theories, we formulate 
the junction conditions which allow for discontinuities in the first derivatives of the metric. 
In Sec. m we consider f{R) theories in the Einstein frame, i.e., Einstein-scalar theories 
conformally equivalent to f{R) theories, and formulate the junction conditions in the Einstein 
frame. In doing so, we suggest a new type of gravitational coupling for the matter on the 
brane. Finally, in Sec. [5l we translate these generalized junction conditions back to the 
original, Jordan, frame. We work in D spacetime dimensions. The gravitational constant 
SttG in D dimensions is set to unity; the signature is — h + • — h. 

§2. "Weak" junction conditions in quadratic theories of gravity 

Consider the quadratic Lagrangian, 

L = -2A + R + -f RabcdR^''^'' - 4/3 RabR^"" + aR\ (2-1) 

The variational derivative of y/—gL with respect to the metric g^B yields, up to a divergence, 
6{y^L) = -^/^(T^^ 5gAB with (see e.g. 25)) 

(^AB = ~^ 

+ 27 Ralnp Rb'-'"'' + 4(2/3 - ^)R''d R^'bac - ^1 Ra"" Rcb + 2a R Rab ^^-2) 
+ 4(7 - P)nRAB + 2(a - l3)gABnR- 2(« + 7 - 2P)DaDb R 
where D is the covariant derivative associated with gAB- The field equations are, in the bulk, 

(TAB = 0. (2-3) 

These equations are fourth order in the derivatives of the metric, except for the Gauss- 
Bonnet combination a = /? = 7 that we shall exclud^. 

We use Gaussian-normal coordinates in which the metric is 

ds^ = dy^ + 7^^dx^da;'' , (2-4) 

where the brane is assumed to be located at ?/ = 0. In terms of the extrinsic curvature, 

AV--i^. (2-5) 



The junction conditions in Einstein-Gauss-Bonnet theory are weh-known. See 26), and, e.g. 27). 
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the Riemann tensor is 

OK, 



m = ^ -\- X E = D K — D K 



where Dp and R^^p^j are the covariant derivative and the Riemann tensor associated with 
the metric 7^,^, with the Greek indices being raised or lowered by the metric 7^1,. 
For convenience, we introduce the tensor 



Keeping in cab only the terms with highest derivatives in y, i.e. the terms proportional to 
dHpi,/dy, one finds (see e.g. 25)) that there are no such terms in ayy and ayp, and that they 
appear in a^j, under the following combination : 

cTp. = 4^[(7 - P)Hp, + {a- P)ip.H] + • • • , (2-8) 

with H = Y'^Hp^. 

Suppose now that there exists a sub-class of metrics jpuiy, x^) solving the bulk equations 
o'AB = 0, whose third order derivatives Hpi, jump across y = 0; that is, such that Hpi, can 
be written as, e.g., Hp^ = hpy^x^) tanh(?//£) with £ — >• 0, in a region of order £ around y = 
( "Z2-symmetric" case). For this sub-class of metrics apy exhibits a Dirac distribution- like 
behavior at y = : 

4^ [(7 - P)Hp, + (a - P)^p,H] = 5{y)Dp, , (2-9) 
where Dp^, is the "strength" of the singularity. Integration across the brane then yields 

4[(7 - (3)Hp, + (« - (3)^p,H]t = Dp, , (2-10) 

where = limj^^o+ Fiy) — hnij^^o- F{y) . If we require this class of metrics to be solutions 
of the field equations, then we must have 

Dp, = Spy, (2-11) 

where we may naturally interpret 5*^,^ as the total energy-momentum tensor of the brane. 
These equations, together with Eq. fl2-10p . give the junction conditions. 

These junction conditions were first given in 28) (in the particular case D = A and hence 
7 = 0, because of the Gauss-Bonnet theorem) and generalized recently in 29) to Lagrangians 
of the type L = f{R,RABR^^,RABCDR'^^'^^)- It is clear that such discontinuities are 
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specific to higher derivative theories and smoothly disappear in the Einstein hmit when the 
parameters a, (3 and 7 are "switched off." 

An important point (apparently overlooked in 29)) is that the contracted Bianchi iden- 
tities, DbCt^a = 0, imply that the brane energy-momentum tensor is conserved: 

D,S\ = 0. (2-12) 



We leave to further work the question of finding a bulk metric, a solution to the bulk 
equations (e.g. anti-de Sitter), which can be written in Gaussian coordinates in such a way 
as to exhibit a discontinuity in its third derivative across y = 0. (The forms of the metric 
given in, e.g., 30) do not belong to this desired sub-class as they exhibit a jump in their first 
derivatives.) 

A few remarks to conclude this section are in order : 

(a) For 7 = 2/9 — a (generalization of Eddington's choice^^^), the junction conditions (12- lip 
take the form (in the Z2-symmetric case) 

(b) For a — f3 = —-j^ (that is Weyl's choice in D = 4 and with 7 = 0^^), they become 

Therefore the total energy-momentum tensor on the brane must be traceless. 

(c) Finally in the pure f{R) case (/? = 7 = 0, L = —2A + R + aR"^), the total energy- 
momentum tensor on a Z2-symmetric brane is constrained to be of the form 

S^, = 8aH^^,. (2-15) 

The fact that S^u is conserved implies if is a constant. Thus, the matter on the brane 
must be vacuum energy. 

This final example indicates that imposing the metric to be of class (continuity of 
7^1. and its first and second derivatives) is probably too restrictive to allow for physically 
interesting braneworld scenarios in higher derivative theories of gravity. We shall therefore 
seek junction conditions that allow for more singular metrics. In the rest of the paper, we 
concentrate on f{R) theories. 
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§3. Junction conditions in the Jordan frame : the standard approach 

The field equations derived from the variation of y/—gf{R) with respect to the metric 
gAB are cab = in the bulk, with 

(TAB = f\R)GAB + \gAB{Rf\R) - fm + QABafiR) - DADsfiR) , (3-1) 



where f'{R) = df /dR and Gab = Rab — ^9abR is the Einstein tensoiCj. In a Gaussian 
normal coordinate system, ds^ = dy"^ + 'y^i,dx^dx'^ , the Einstein tensor is decomposed as 

Gyy = -^{K^^K^'' ~ + R) , 

Gy, = -D,{K\ - 6;K) , (3-2) 
G^, = dy{K,, - K-f^,) + 2K%Kp, - 3KK^, + ]p^,{K,,^K^^ + K^) + G,, , 

where the extrinsic curvature is K^^ = —\dy'y^y. We note that 

R = 2dyK - K^^K^"" -K^ + R 

D (3-3) 

= 2dyK - k;,k*^-^ - -^—^K^ + R , 

where K*^ is the traceless part of the extrinsic curvature : K*^^ = K^^ — j^j^u- 
For convenience, we decompose aAB into 

= Qab + Lab , (3-4) 

with 

Qab = f\R) Gab + \{R f\R) - f{R)) gAB , 

2 (3-5) 

Lab = -DADsfm + gABnf\R) . 
Their components are given by 

Q^^ = ^(Rf(R)-f(R))+f'(R)Gyy, 

Qy, = f'{R)Gy, , (3-6) 
Q,u = \i,u{Rf\R) - fm + f\R)G,, , 

Lyy = -Kdyf\R)+nf{R), 

Ly, = -d,dyf\R) - K;dJ\R) , (3-7) 

V = -D,.f{R) + K,,dyf\R) + l^,{dyyf\R) + Uf'{R) " K dy f' {R)) . 



and 



We consider here metric f{R) theories. For variation a la Palatini, see 32) and references therein. 
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In Sec. |2]we considered the class of metrics 7^,^ which were continuous across y = with 
continuous first and second derivatives. In that case all the components of the Einstein tensor 
and of Qab are well behaved and a Dirac distribution appears in L^j, in the term '~i fiudyy f {R) . 
Now, since, for that sub-class, R is continuous and hence dyR is at most discontinuous, the 
delta distribution behavior of dyyf'{R) comes from dyyR: 

dyyf'iR) = f"{R)dyyR + rmdyRf 

= 2f"{R)dyyyK + " ■ " = 2f"{R)dyH + > > > , 

with = dyyK^y. In the particular case f{R) = —2A + R + aR'^ we thus recover the result 

of Sec. m 

Here, on the other hand, we shall consider the class of metrics which are continuous 
across y = but which allow for (a certain type of) discontinuity in their first derivatives. 

The scalar curvature R which could be now, a priori, a delta function, must be at most 
discontinuous; otherwise unacceptable terms would appear in Qab (unless of course 

f{R) = —2A + R). Furthermore, an inspection of L^^ tells us that R must be in fact 
continuous across y = 0; otherwise a term would appear in L^i, (unless f{R) is 

quadratic, because it precisely arises from the term f"'{R){dyRY , see (13-81) ). We shall 
therefore restrict ourselves to the sub-class of metrics with continuous scalar curvature R. 

Again we leave to further work the question of finding a bulk metric, solution to the 
bulk equations, which can be written in Gaussian coordinates in such a way as to exhibit 
discontinuities in its first order derivative across y = 0, but not in R. (It seems that the fact 
that the scalar curvature must be continuous has been overlooked in the recent paper. ^^^) 

We then see by inspection that the yy and yfi components of ctab at most jump across 
y = and that the delta-like part of a^u is 

dy[f'{R){K^, - K^,,) + -i,J"{R)dyR] = 5{y)D,, . (3-9) 

Integration across the brane then yields the junction conditions 

D^, = [nR){K^, - K-f^,) + ^,J"{R)dyR]t = S,, , (3-10) 

where S^j^y is the total energy momentum tensor on the brane. From the contracted Bianchi 
identities, we have that it must be conserved: D^S'^fj^ = 0. 

Note that when f{R) —>■ R the junction conditions (13-101) do not reduce to the familiar 
Israel conditions'^-* as they have to be supplemented by the condition of continuity of R. 
What happens when f{R) = —2A + R + i'^R'^ + ... when ^ ^ is that the bulk geometry may 
approach a solution of the Einstein bulk equations (e.g., AdS) everywhere, to the exception 
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of a region of size £ in the vicinity of the brane, so that when £ becomes very small the thin 
shell limit is no longer valid and the thickness of the brane must be taken into account. To 
render this irreducible difference between Einstein and f{R) theories manifest, let us split 
the junction conditions into their trace and traceless parts, recalling that 7^,^ as well as R 
are continuous (see (13-31) ) : 



[R]t = =^ [K]t = 0, [2dyK - K;^K*n^_ = , ^^^^^^ 



f'{R)[K;^]t = s;^, 
iD-i)riR)[dyR]1: = s, 

where K*^, and S*,^ are the traceless parts of the extrinsic curvature and brane energy- 
momentum tensor, respectively. The "weak" junction conditions considered in Sec. [2] are 
just a particular sub-class of the above such that 5**^, = 0. 

For further reference, we shall also generalize them to a non-Gaussian coordinate system, 
ds^ = N'^dy'^ + 'y^^dx^dx'^ , where is a continuous lapse function: 

Mt = 0, [K]t = 0, [R]t = 0, 

f'{R)[K;^]t = s;^, (3.12) 

{D-l)nR)j^[dyR]t = S, 

where, now, the extrinsic curvature is defined as Kf^i, = —jf^dy'-y^u- 

In this and the preceding section we just stated that the jumps in some derivatives of 
the bulk metric coefficients must be equal to the total energy-momentum tensor of matter 
on the brane. We have yet to decide on the type of matter we want to have on the brane. 
An even more crucial issue is to decide how matter on the brane couples to gravity, that is, 
how it couples to the metric 7^,^ (which is the only dynamical variable at our disposal). 

It is natural (and this is the hypothesis which is "universally" made in the literature) to 
assume that matter is minimally coupled to the metric. For example, if matter on the brane 
is taken to be a scalar field with potential V{iIj), we have 

S,. = d.^d^^P - 7^. (^^r'^d^ijdp^ + V{^)^ . (3-13) 

As we shall see in an accompanying paper, '^^^ the junction conditions (13-111) and (I3-I2p thus 
obtained are the standard ones, that is, those one derives from the action ^ J d^x \/—gf{R) 
supplemented by the Hawking-Luttrell boundary term J d^^^ Xy/ —'j f (R) K .^^^ 

In the next section, we shall question this conventional wisdom, taking advantage of the 
conformal equivalence of f{R) theories of gravity with General Relativity coupled to a scalar 
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field. This will lead us to treat the scalar curvature as an independent field and allow us to 
propose more general junction conditions than (13- lip and fl3-12p . 



§4. Junction conditions in the Einstein frame : extra degree of freedom and its 

coupling to matter on the brane 

In order to simplify the notation, we denote quantities in the Jordan frame with tildes, 
e.g., R for the scalar curvature, and those in the Einstein frame without tildes. Thus all the 
quantities that appeared in the previous section should be tilded. 

As is well-known, the bulk f{R) field equations for the Jordan frame metric (Jab, 
i.e., aAB = with aAs given by (13- 41) . are equivalent to bulk Einstein equations for the 
"Einstein frame" metric, 

gAB = Qab exp I , '^^ I , (4-1) 

with a scalar field minimally coupled to gravity : 

1 



Gab = dA<pdB<p - Qab [^^dc 09^0 + W{<P)j , (4-2) 
where the potential W{(f)) is implicitly defined as a function of cf) via 



W{R) = \{Rf\R)-f{R))f\R)--'\ 0=^§3^1nm. (4-3) 

Because of the Bianchi identities, the Einstein equations (14- 21) are consistent only if cj) satisfies 
the Klein-Gordon equation, 

□^-- = 0, (4.4) 

Mathematically, this conformal transformation transforms the original, fourth-order dif- 
ferential equation (1341) into two second-order differential equations, one for the Einstein 
frame metric Qab-, the other for 0. Physically, it shows that f{R) is a "scalar-tensor" theory 
of gravity where 0, that is, the bulk Jordan frame scalar curvature i?, is an extra, indepen- 
dent, degree of freedom.'^^^ What we shall dwell upon in the following is the coupling of this 
extra degree of freedom with matter on the brane. 

Again we shall use Gaussian coordinates, that is 

= dz^ + 7^^dx^dx^ , (4-5) 

and consider a class of metrics which are continuous across z = 0, but which allow for 
discontinuities in their first derivatives. We shall also impose the scalar field to be continuous, 
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and allow for a discontinuity in its first derivative. Now, is directly related to the scalar 
curvature of the bulk Jordan frame metric ; one must not however deduce hastily that the 
condition of continuity of (j) is equivalent to imposing the continuity of the scalar curvature 
of the Jordan frame metric ; indeed the relation between and R holds in the hulk only and 
they may differ, as we shall see, by a term confined on the brane. 

Thus, allowing for discontinuities of the 2;-derivatives of the metric and of 0, the right- 
hand side of the field equations jumps at most. The G^z and Gz^ components of Gab also 
jump at most. As for the delta-like part of the Einstein tensor G^jy, it is, see (13 -2^ : 

dz{K^,-K^^,) = 6{z)D^,. (4-6) 

Integration across the brane then yields the Israel junction conditions^^^ 

D^, = [K^,-K^^,]t = T^u. (4-7) 

where T^,y is the total energy-momentum tensor of the brane in the Einstein frame. Later 
we shall relate it to the energy-momentum tensor of the brane in the Jordan frame S^^. 

Since the first derivative of is allowed to be discontinuous, 00 also exhibits a delta 
function-like behavior. From the Bianchi identities, 

= dB<P ( 00 - ^ J + Da{T^b6{z)) , (4-8) 

we have 

d^v[dzCt>\t = ~D,T\, (4-9) 

where ip{x^) = (j){z = 0,x^). 

Just as in the case of working in the Jordan frame, the last task is to express the total 
stress-energy tensor of the brane matter in terms of the matter variables. Since the gravita- 
tional variables are 7^jy and (p we have to decide how matter couples to those gravitational 
fields. For example, if matter on the brane is taken to be a scalar field ip with potential 
V{iIj), we may consider as the matter action, 

Sm = -j d^-ia;v^ (^F,{^)^r''d,ijd,tP + F2{^)V{ij)^ , (4-10) 

where Fi and F2 are two a priori arbitrary functions of (p. The associated energy-momentum 
tensor is 

T^, = F^{^)^^iJ^,^P - 7^. Qfi(<^)7"^9«7/>9;3^ + F2(^)l^(^)^ . (4-11) 
The field equation for ip is 

D''{Fi{ip)d,^) - F2{v)V'{^) = . (4-12) 



11 



Thus the divergence of the energy-momentum tensor gives 

D.T\ = -d,^ (^^^"""^"^^/^^ + ^^(^)) • ^4-^3) 

Following Einstein's suggestion to Nordstrom, we require that the source for (j) in (14- 9p 
be related to the trace of the matter energy-momentum tensor, 

D,T\^d^^T. (4-14) 

This imposes 

Fi(<^) = exp[(Z} - ?,)k{^)] , F2(<^) = exp[(D - l)k{^)] , (4-15) 

where k{Lp) still has to be determined. Plugging these back into fl4-10l) . we see the meaning 
of the condition (14- 141) . Namely, the matter should be minimally coupled to a metric 7^,^, 

^m = 5m[V;^], V = e''^''V> (4-16) 

where ip now represents general matter variables not restricted to a scalar field. 
The junction condition fl4-9p then becomes 

(\k 

[dAt = -^T, (4-17) 



with 



T = -2- 



-7 57 



(4-18) 



where the second line is the case when the matter is a scalar field. The junction condition 
fl4-17l) is nothing but the one used when studying the brane cosmology with a bulk scalar 
field.39) 

We also note that the matter action can then be rewritten in terms of the Jordan metric 

Ifxv as 

with 

C{ip) = kiip) + = . (4-20) 

This is where the fact that we are treating a f{R) theory in the Einstein frame comes into 
play. Indeed, if we impose the matter on the brane to be minimally coupled to the Jordan 
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metric, then we must choose (a point aheady known to Einstein^^-* (see also 41))) C = 0, 
that is, 

k{ip) = . (4-21) 

VP-l)P-2) 

This reduces the junction conditions to 

We shall now translate back to the Jordan frame the generalized junction conditions (14- 7p 
and (14- 171) . and show that they reduce to those obtained in f l3-12p with f l3-13p only when 
k{ip) is imposed to be given by fl4-2ip . 

§5. Back to the Jordan frame : generalized junction conditions 

In a Gaussian normal coordinate system in which the line element reads ds^ = dz^ + 
7^!ydx'^da;'^, the junction conditions we have obtained in the Einstein frame are 

Ml = o, m; = o, 

where K^^ = —^d^'y^u and the second line recalls that the induced metric as well as have to 
be continuous across the brane. We assume that the matter on the brane couples minimally 
to the metric ^y^^ = e^'^^'^^jf^i^ as given by fl4-16p . 
Let us perform the following transformations 



^In/'(P) 



W{4>) ^ w(p) = -(p/'(p) - f(p))f(p)-T^' , f^'^) 



2 

-2 



gAB gAB = gABf'ip) 

It is a side exercise to show that if gAB{,x^) and 0(a;'") are solution of the bulk field equations 
(14- 2p then p{x^) is the scalar curvature R of the bulk Jordan metric cjAB- However if one 
includes the presence of matter on the brane, there appears a delta function-like singularity 
in R while p is continuous unless the matter is minimally coupled on the brane, i.e., unless 
C{ip) = (modulo a constant), as we shall see below. 

Note that the coordinates are no longer Gaussian for the Jordan line element 

ds' = gABdz^dx"" = f'{p)-T^^gABdx^dx'' = f'{p)-T^^dz^ + ^.dx-dx" , (5-3) 
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with 7^1, = /'(p) ^'^'Jfj.u the induced metric of the Jordan brane. We introduce 

Kf.u = -lf'^'d,%,, k = r''K,,- (5-4) 

The junction conditions (15-11) then translate as follows. As already mentioned, the con- 
tinuity of = f'(p) translates into the continuity of p and the continuity of the 
metric qab translates into the continuity of the Jordan induced metric : 

[p]l = 0, [%ut = 0. (5-5) 

As for the jumps in and the extrinsic curvature they translate into 

where 

~ 2 ^g^[e^^(^)7p.;^] 

and a star means taking the traceless part. We note that because of the assumption that 
the matter is minimally coupled to the metric 'j^^, we have 

yZ^T = -2^" e-- = ^/^T. (5-8) 

If matter on the brane is a scalar field ip, then is given by fl4-17l) so that 



T = e^ 



(D-i)c{^) fe~2C(^)^^^p^,^p - f^e-'^^i^^^'^Op^d^^ + V^] ] . (5-9) 



Now, generalizing fl3-3p to a non-Gaussian coordinate, the scalar curvature R may be 
expressed as 

R = 2N-^d,k - kp^k^''' -k^ + R; N = f\p)~T^ . (5-10) 
Integrating this across z = Q and using the junction conditions (15-61) . one finds that 



' n _ 1 fir* T 

fliVd. = 2|A1! = -2,/— --. (5.11) 



From this, we deduce that 



/ n _ 1 An 1 T 
R = p-2J- — -^^fo^^Ls^z). (5-12) 
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Therefore, as anticipated, the continuity of translates in the continuity of p but not of R . 

The junction conditions fl5-5l) and fl5-6p are the central result of this paper. When the 
arbitrary function C vanishes, we have from fl5-12p that R = p everywhere including on the 
brane and we have that T^,y coalesces with 5*^1,, the stress-energy tensor of matter minimally 
coupled to the brane metric introduced in f l3-10p . The junction conditions thus reduce to 
those obtained in Sec. [31 

When C is a non-trivial function, they generalize them to the case when matter on the 
brane is coupled not only to the brane metric 7^1, but to the extra degree of freedom of f{R) 
gravity which, in the Jordan frame, is the quantity p, equal everywhere but on the brane to 
the scalar curvature R. 

We just note here, to conclude, that, if we choose C{(p) = (p/^/{D — 1){D — 2), that is 
k{(p) = 0, then the junction conditions in the Jordan frame closely resemble the standard 
Israel junction conditions 

f[K,,,-K^,,,]t=f^.. (5-13) 
However the coupling of matter to the brane metric is not minimal, as given by {\5-7\i . 

§6. Conclusion 

We thoroughly investigated the junction conditions in f{R) theories of gravity. We found 
that in a pure f{R) theory in which matter on the brane couples minimally to the metric, the 
bulk scalar curvature R must be continuous across the brane, which is in marked contrast 
with the case of Einstein gravity. Then taking advantage of the conformal equivalence of f{R) 
theories with Einstein gravity with a scalar field, we clarified the importance of identifying 
the scalar curvature R as an extra degree of freedom and the specific form of the coupling 
of matter to this extra gravitational degree of freedom on the brane. 

Then as a bonus of working in the Einstein frame, we presented a natural generalization 
of the coupling of matter to gravity. In the original frame, this leads to a non-trivial coupling 
of the matter on the brane to the extra degree of freedom, which allows a delta function-like 
behavior of the scalar curvature. This suggests a new class of braneworld models whose 
solutions may have a smooth limit in the Einstein limit f{R) — > R- 

It is known that, in the bulk, an f{R) theory may be rewritten as a Brans-Dicke theory 
with uj = but with a potential. If we use this equivalence, the generalization mentioned 
above may be regarded as a non-trivial coupling of the Brans-Dicke scalar to the matter on 
the brane. 

We leave to an accompanying paper the derivation of the junction conditions (15-60 via 
a first order description of the f{R) action as well as an analysis of the braneworld models 
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that they may lead to. 
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